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Abstract 

Heterozygote disadvantage is potentially a potent driver of popu¬ 
lation genetic divergence. Also referred to as underdominance, this 
phenomena describes a situation where a genetic heterozygote has a 
lower overall fitness than either homozygote. Attention so far has 
mostly been given to underdominance within a single population and 
the maintenance of genetic differences between two populations ex¬ 
changing migrants. Here we explore the dynamics of an underdom¬ 
inant system in a network of multiple discrete, yet interconnected, 
populations. Stability of genetic differences in response to increases 
in migration in various topological networks is assessed. The network 
topology can have a dominant and occasionally non-intuitive influence 
on the genetic stability of the system. 

Keywords; Under dominance. Coordination Game, Network Topology, 
Dynamical System, Population Genetics 


1 Introduction 

Variation in the fitness of genotypes resulting from combinations of two al¬ 
leles {e.g., A- and B-type alleles combined into AA-, AB-, or BB-genotypes 
resulting in waa, wab-, and wbb fitnesses respectively) result in different 
evolutionary dynamics. The case in which a heterozygote has a lower fitness 
than either homozygote, wab < WAA and wab < wbb-, is termed under¬ 
dominance or heterozygote disadvantage. In this case there is an internal 
unstable equilibrium so that the fixation or loss of an allele depends on its 
starting frequency. In a single population, stable polymorphism is not possi¬ 
ble. However, when certain conditions are met, populations that are coupled 
by migration (the exchange of some fraction of alleles each generation) can 
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result in a stable selection-migration equilibrium. This selection-migration 
equilibrium is associated with a critical migration rate [m *); above this point 
the mixing between populations is sufficiently high for the system to behave 
as a single population and all internal stability is lost (Altrock et al. 2010). 

Under dominance can be thought of as an evolutionary bistable switch. 
From the perspective of game-theory dynamics it can be interpreted as a 
coordination game (Traulsen and Reed, 2012). The properties of under¬ 


dominance in single and multiple populations have led to proposals of a 
role of underdominance in producing barriers to gene flow during speciation 


(Faria and Navarro, 2010; Harewood et ah, 2010) as well as proposals to 


utilize underdominace both to transform the properties of wild populations 


in genetic pest management applications (Curtis 

1968 

Davis et al. 

2001 

Sinkins and Gould, 20061 Reeves et al. 

2014 

) and to engineer barriers to gene 


flow (transgene mitigation) from genetically modified crops to unmodified 


relatives (Soboleva et ah, 2003 Reeves and Reed, 2014) 


The properties of underdominance in a single population are well under¬ 


stood (Fisher, 1922; Haldane, 1927; Wright, 1931) and the two-population 


case has been studied in some detail (jKarlin and McGregor| 

1972a||b 

Lande 

1985 

Wilson and Turellil 1986| Spirito et al. 19911 Altrock et a 

12010 


2011 


with fewer analytic treatments of three or more populations (Karlin 


and McGregor 1972a b). Simulation-based studies have been conducted for 


populations connected in a 

lattice ( 

Schierup and Christiansen 

1996 

Paynej 

et al. 

2007 

Eppstein et al. 

2009 

Landguth et al. 

2015 

) and “wave” approx- 


imations have been used to study the flow of underdominant alleles under 


conditions of isolation by distance (Fisher, 1937 Pialek and Barton, 1997 


Soboleva et ah, 2003 Barton and Turelli, 2011). Despite these properties 


and potential applications, underdominance has been relatively neglected 
in population genetic research (Bengtsson and Bodmer 1976). A large fo¬ 


cus of earlier theoretical work with underdominance was on how new rare 
mutations resulting in underdominance might become established in a pop- 


ulation ( 

Wright 

TMTl 

Bengtsson and Bodmer 1976 

Hedrick 

TMTj 

Walsh 

19821 Hedrick and Levin, 119841 Lande, 1984, 19851 

Barton and Rouhani 

19911 Spirito 

1992 

). However, here we are addressing the properties of how 


underdominant polymorphisms may persist once established within a set of 
populations rather than how they were established in the first place. 

We explicitly focus on discrete populations that are connected by mi¬ 
gration in a population network. We have found that the topology of the 
network has a profound influence on the stability of underdominant poly¬ 
morphisms that has been otherwise overlooked. This influence is not always 
intuitive a priori. These results have implications for the effects of network 
topology on a dynamic system (see for review Strogatz[ 2001), particnlarly 
for interactions related to the coordination game (such as the stag hunt, 
Skyrms, 2001), theories of speciation, the maintenance of biological diver¬ 


sity, and applications of underdominance to both protect wild populations 
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from genetic modification or to genetically engineer the properties of wild 
populations—depending on the goals of the application. 


Methods and Results 


We are considering simple graphs in the sense of graph theory to repre¬ 
sent the population network: each pair of nodes can be connected by at 
most a single undirected edge. A graph G = is constructed from 

a set of nodes, ^ (also referred to a vetexes), and a set of edges, S’, that 
connect pairs of nodes. For convenience V = |o/F| and E = \S\, we chose 
V (for vertex) to represent the number of nodes to avoid future conflict 
with N symbolizing finite population size in population genetics. A node 
corresponds to a population made up of a large number of random-mating 
(well mixed) individuals (a Wright-Fisher population, (]Fisher 1922, Wright 


1931) with independent Hardy-Weinberg allelic associations, (Hardy, 1908)) 


and the edges represent corridors of restricted migration between the pop¬ 
ulations. We are also only considering undirected graphs: in the present 
context this represents equal bidirectional migration between the popula¬ 
tion nodes. Furthermore, we are only considering connected graphs (each 
node can ultimately be visited from every other node) and unlabeled graphs 
so that isomorphic graphs are considered equivalent. 

The network graph G is represented by a symmetric V x V adjacency 
matrix £/. 



« 1,1 

0 - 1,2 

« 1,3 

■ ■ ■ Oiy 


« 2,1 

02,2 

« 2,3 

■ ■ ■ 02,V 

= 

« 3,1 

03,2 

03,3 

■ ■ ■ 02,V 


av,i 

OV,2 

OV,3 

. . . ov,v_ 


The presence of an edge between two nodes is represented by a one and 
the absence of an edge is a zero. The connectivity of a node is 



j&.yr 


Each generation, g, the allele frequency, p of each population node, i, is 
updated with the fraction of immigrants from n adjacent populations, j, at 
a migration rate of m. 


N 

Pi^g = (1 - Cim)pi^g-i -f 

i=i 
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Note that this equation will not be appropriate if the fraction of alleles 
introduced into a population exceeded unity. See the discussion of the star 
topologies illustrated in Figure 

The frequencies, adjusted for migration, are then paired into genotypes 
and undergo the effects of selection. The remaining allelic transmission sum 
is normalized by the total transmission of all alleles to the next generation 
to render an allele frequency from zero to one. 

, ^ pjg + - Pi,9) 

pig + - Pi,g) + (1 - Pi,g)'^ 

Note that here for simplicity we set the relative fitness of the homozy¬ 
gotes to one, waa = = 1 and the heterozygote htness is represented 

by wab = w. 


2.1 Analytic Results 


Under dominance in a single population has one central unstable equilibrium 
and two trivial stable equilibria at p = 0 and p = 1. When one considers 
multiplying the three fixed points of a single population into multiple di¬ 
mensions it can be seen that, if migration rates are sufficiently small, nine 
equilibria (3x3) exist in the two-population system. Again, there is a central 
unstable equilibrium, the two trivial stable equilibria, and two additional in¬ 
ternal stable equilibria. The remaining four fixed points are unstable saddle 


points that separate the basins of attraction (see Figure 3 of Altrock et al. 


2010). However, as the migration rates increase the internal stable points 


merge with the neighboring saddle points and become unstable themselves 
(Karlin and McGregor, 1972b Altrock et ah, 2010). In three populations 


there are a maximum of 27 (3 x 3 x 3) fixed points. The three types found 
in a single population, six internal stable equilibria, and 18 saddle points. 
The general pattern is that in V populations there are a maximum of 3'^ 
equilibria if migration rates are sufficiently low. There will always be one 
central unstable equilibrium (a trajectory starting near this point can end 
up in any of the basins of attraction) and two trivial fixation or loss points. 
A U-dimensional hypercube represents the state space of joint allele frequen¬ 
cies of a U-population system. This hypercube has 2^ vertexes (corners). 
It can be seen that the internal stable points, if they exist, are near these 
corners when the corners represent a mix of zero and one allele frequencies. 
There are exactly two corners, the trivial global fixation and loss points, 
that do not contain a mix of unequal coordinate frequencies. Thus, there 
are 2^ — 2 possible internal stable points that represent alternative conhg- 
urations of a migration-selection equilibrium polymorphism. Finally there 
are up to 3'^ — 2^^ — 3 saddle points that separate 2^ possible (but never 
less than two) domains of attraction. 
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We have a system of equations, p[ g, that describe the dynamics of allele 
frequency change within each population and these dynamics are coupled by 
migration. It is useful to think of the difference in frequency each generation, 
^P,i — Pig ~ Pi g-i- We can set = 0 to solve for fixed points in the state 
space. However, the system needs to be simplified in order to be tractable. 
For example if we only look along the pi = p 2 axis in the two-population 
case we get three solutions, p = 0, p = 1, and 


PWie-yy = 


waa - wab 


waa + wbb - ‘^wab 


in the general case and 


Pii^.jV 


1 — U) 
2 - 2a; 


in the simplified (equal homozygote fitness) case. The first two points are the 
trivial loss of polymorphism. The third is identical to the internal unstable 


equilibrium point in a single population (Altrock et ah, 2010). 


In fact this unstable equilibrium point is always found along the pi = 
P2 = ■ ■ ■ = Pv axis. Note that the migration rate, m, is not a part of the 
solution. The position of this point in the state space is independent of 
migration rates. Since it falls along the axis where the allele frequencies of 
all populations are equal, migration between populations, and in fact the 
population network topology itself, has no effect. 

While the position of this internal point remains fixed regardless of the 
number of interacting populations in the network, there is a multiple pop¬ 
ulation effect on the rate of change away from this point. Solving for the 
eigenvalues, Aj, of a Jacobian matrix, J, of partial derivatives of the system 
along the Pi = P2 = ''' = Pv axis shows that the rate of change follows the 
pattern 

21 / 


^ilp=i - 


-V. 


1 -h UJ 

Thus, as the number of interconnected populations increase the magnitude 
of flow away from the internal point increases. This rate is a function of 
both the heterozyote fitness (a;) and the number of interlinked populations. 
At a; = 0, or lethality of the heterozygous condition, the rate of change is 
equal to the number of connected populations. 

We are also interested in the internal stable equilibria that allow dif¬ 
ferences in allele frequencies among populations to be maintained. In the 
two-population special symmetric case this can be solved from J along the 
Pi = 1 — P2 axis and yields 


A2 — — 


1 -Lu + 2m*(2m* -3) 

(1 — 2m*y 


m 


A 2 


=0 = - (3 - x/5T4J) 
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Figure 1: The stability of some simple network configurations. For each 
number of nodes and network topology (linear in blue, cyclic in red, star- 
like in yellow, and fully interconnected in green) the critical migration rate 
allowing polymorphic underdominant polymorphism to persist at w = 1/2 
is plotted. Examples of each type of network at E = 5 are plotted as 
graphs in the legend to the right. The shading of the nodes represents the 
allele frequency between zero and one of each population near the critical 
migration rate value. 


(see Appendix A of |Altrock et al. 2010, for more detail). Unfortunately, with 
three or more populations, even with highly symmetrical configurations, we 
have not found a single axis or plane through the state space that captures 
these internal stable points. Therefore, we have used numerical methods to 
characterize the critical points. 


2.2 Numerical Simulations 

Sets of populations were initialized with approximately half (depending on 
if there were an even or odd number of populations) of the nodes in the 
network at allele frequencies near zero and half near a frequency of one. The 
allele frequencies were offset by a small random amount from zero or one 
to avoid being symmetrically balanced on unstable trajectories. Migration 
rates were slowly incremented at each step and the system was allowed to 
proceed to near equilibrium, when the difference in allele frequencies between 
generations was less than 10“^^, before the next step. The process was 
repeated until the point where a collapse in differences of allele frequencies 
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between all populations was detected. This point was then reported as the 
critical migration rate of the network. 


2.2.1 Example Topologies 


The stability of a range of basic network topologies were investigated, shown 
in Figure In general, in these examples, the diameter of the network is 
predictive of migration-selection stability. Linear configurations had the 
highest stability, cyclic configurations approximately reduced the both the 
diameter and stability in half. Fully interconnected populations had both 
the smallest diameter and lowest stability. 

Another pattern that became apparent is an even-odd alternation in 
stability. Except for starlike networks, an odd number of nodes results in a 
relatively lower stability than an even number of nodes. 

Starlike networks showed an interesting pattern. They were approxi¬ 
mately of the same stability as cyclic networks; however, the even-odd alter¬ 
nation was inverted—odd V graphs showed enhanced stability, showing that 
the even-odd pattern is not absolute. At this heterozygote fitness (cj = 1/2) 
starlike networks with greater than eight nodes could not be evaluated. At 
F > 9 before the critical migration rate is reached the total amount of 
immigration into the central population exceeds 100%. 

Fully interconnected networks had the lowest stability and by far the 
greatest number of edges. Unlike the other networks the fully-interconnected 
systems declined at higher V. However, the number of edges grew much 
faster than the number of nodes. Note that the even-odd alternation in 
relative stability is still apparent, even in these graphs. 

The effects of a range of topologies for six nodes and five edges was also 
explored, shown in Figure All of these networks have identical treeness, 
<^tree = 1 {sensu pCie and Levinson 2007). In general the stability is cor¬ 
related with the diameter of the network. However, the clear exception is 
the “double-Y” topology, which has the highest stability of all. This has 
inspired an alternative measure of the treeness of a network that we will 
refer to as “dendricity” to avoid conflicting with prior definitions of treeness 
in the literature. 
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Figure 2: The stability of all possible simple connected networks made of 
three populations with five corridors of migration. Here the critical migra¬ 
tion rate was evaluated at a relative heterozygote fitness of cu = 3/4. The 
networks are arranged by diameter of the network declining from five on the 
left to two on the right. The shading of the nodes represents the allele fre¬ 
quency between zero and one of each population near the critical migration 
rate value. The shading of the bars is set to 50% gray to aid visualization 
of the allele frequencies of the nodes. 


2.2.2 Random Graphs 

In order to evaluate general correlations between migration-selection stabil¬ 
ity and network summary statistics we generated 100,000 random connected 
graphs of up to 20 nodes in size and evaluated their stability. The results 
are summarised in Supplementary Table 1. We found that the most stable 
network configurations contained nodes with at most three edges (see the 
next section below) so, in order to avoid the problem of the rate of immigra¬ 
tion more than replacing a local population we set a maximum migration 
rate of m = 1/3 and reported this as m* for the subset of highly stable 
network topologies—only 0.64% of the random networks reached this point 
of m = 1/3—these highly stable networks are explored in the next section. 

The following parameters were estimated from these networks: Variance, 
which refers to the variance in connectivity {ci the number of edges incident 
with node i) of all the nodes in the network. Efficiency refers to the shortest 
path lengths between nodes in the network according to 

f^efficiency = ^ 

where d{i,j) is the minimum path length between nodes i and j. The 
diameter of the graph is the maximum d{i, j) G G. Dendricity is the fraction 
of nodes incident with three edges where at least one edge is a “bridge” edge 
(removing bridge edges results in an unconnected graph) out of the total 
number of internal nodes. Evenness is simply a binary variable of zero or 
one to indicate if an odd or even number of nodes are present in the graph. 
Finally, terminalness indicates the fraction of nodes in the graph that are 
terminal (or leaf) nodes. 

Each of the summary statistics we addressed were significant predictors 
of network stability; however, because of correlations between these mea¬ 
sures caution must be used to interpret the results. For example, contrary 
to intuition the number of nodes was negatively correlated with stability. 
This is because the number of possible edges, which generally lower stabil¬ 
ity, increases dramatically with the number of nodes. When the number of 
edges is controlled the number of nodes becomes strongly positively corre¬ 
lated. Simply the number of nodes per edge {V/E) is a powerful predictor of 
stability. In general diameter is a strong predictor of stability, particularly 
if the number of nodes are held constant (compare to Figure]^. Dendric¬ 
ity and terminalness also performed well as general predictors of stability. 
Evenness continues to be a predictor of stability, specifically even ordered 
networks have a greater stability than odd ordered ones, but the predic¬ 
tive power is weak compared to other measures. Variance and efficiency 
are a bit more difficult to understand. Increased variance in the number 
of edges per node is associated with lower stability, while intuitively one 
might expect the opposite. When this is measured as variance divided by 
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the total number of edges the correlation almost disappears and in fact be¬ 
comes slightly positive—as expected—unless the number of nodes are held 
constant. Efficiency also varies in a non-intuitive way. It is positively cor¬ 
related unless either the number of nodes or the number of edges are held 
constant where it becomes strongly negatively correlated; however, if both 
are constant efficiency becomes slightly positively correlated again. 


When exploring model selection to find the minimal adequate model to 
predict m* via adjusted and Mallows’ Cp as implemented in the R pack¬ 
age “leaps” all ten summary statistics were retained using all four meth- 


ods (R Development Core Team 

2008 

T. Lumley using Fortran code by 

A. Miller 

2009 

). Using all of the predictors in the full linear model explains 


the majority, = 0.72, of the variation in m*. 


2.2.3 Evolving Networks 

In order to more fully explore the upper edge of highly stable networks for 
a given V we wrote a program that would make random changes to the 
network and evolve higher stability configurations. Starting from a fully 
interconnected network with as close to half of the nodes near an allele 
frequency of zero or one as possible, edges were randomly selected to be re¬ 
moved or added with the constraint that the new network remain connected. 
Most frequently a single edge was altered but with reducing frequency two 
or more edges could be changed simultaneously to allow larger jumps in 
topology and movement away from locally stable configurations. When a 
network is altered its value is determined. If the new critical value is 

higher than the value of the current network the new network is 

adopted for the next step. If the new critical value is lower than the current 
network, the new network is adopted with a probability equal to the ratio 
of the new and current critical migration values {mnew/''^current)- This also 
allows the network evolution to explore regions off of local maxima. Up to 
V = 5 the most stable network configuration was a linear topology. From 
U > 6 networks with greater stability than the linear configuration were 
found and are illustrated in Figure Note that the number of possible 
connected networks increases dramatically with larger V. The most stable 
networks found in Figure are not expected to result from an exhaustive 
search, particularly for V > 10. They are however meant to illustrate some 
general properties of highly stable networks. 

2.3 Software Availability 

All simulations for both random and evolving networks were written in 
Python 2.7.10. The code is freely available on GitHub: https://github. 
com/akij arl/NetworkEvolve 
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Figure 3: A plot of the stability of highly evolved networks (w = 0.95). For 
comparison the stability of corresponding linear structures is also plotted 
(linear, blue). The most stable network found for each V is given (at the 
end of 10,000 steps of random changes with a Metropolis-Hastings-like chain 
update, starting from a fully interconnected network, over 8 independent 
replicate runs), even to the left and odd to the right near each corresponding 
plotting point. With small change, networks can be substantially more 
stable than the linear configuration. This seems to derive from a balance of 
increasing diameter and forked anchoring structures at either end. 
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3 Discussion 


One result that is beginning to emerge from the study of evolutionary dy¬ 
namics on graphs is that the resulting properties can be sensitive to the 


network topology, but often in a non-intuitive way, {e.g. Hindersin and 


Traulsen, 2014, 2015). There are some general factors that influence, or 


are predictive, of the stability of underdominant polymorphisms in a popu¬ 
lation network. In smaller networks the influence of diameter and evenness 
are apparent. The larger the diameter of the network the effectively lower 
the migration rate is between the edges of the system, because alleles have 
to be exchanged via intermediate nodes; it is well understood that a lower 
migration rate enhances migration-selection stability (Altrock et al. 2010). 
In contrast the more edges there are in a system the higher the effective 
migration rate across the network, which results in lower stability. 

The role of the even-odd number is nodes is more subtle. In a system 
of coupled populations having an allele frequency near p = 1/2 is inherently 
unstable with underdominant htness effects. An odd ordered network that 
is anchored at high and low frequencies near its edges pushes the central 
population near p = 1/2, which destabilizes the entire system. Of course 
there are exceptions to this rule such as the star network topology (star 
networks also have other unusual properties such as acting as amplifiers of 


selection Frean et ah, 2013; Adlam et al. 2015 Hindersin and Traulsen 


2015), and the importance of evenness declines with larger networks. This 


pattern may also change if there are three or more alleles that are underdom¬ 
inant with respect to each other. Interestingly, this effect seems to have been 
completely overlooked in previous work using either numerical techniques or 
wave approximations to study underdominant-like effects. 

The networks that were evolved to higher stability illustrate the effects of 
having two “anchor” nodes at each end of a central linear network “trunk.” 
The anchors are made up of variations on a theme of a node with three edges 
with at least one of these edges being a bridge edge. Intuitively, the flow of 
one allele along two paths into a population can overcome the flow of the 
alternative allele along a single path, thereby enhancing the stability of the 
anchored edges of the system (c/. the discussion of “stem” structures, reser¬ 
voirs, and the movement of clusters of mutants within “superstar” networks 
in Jamieson-Lane and Hauert, 2015). In contrast a strictly linear network 


allows adjacent populations to collapse one by one without any local restric¬ 
tions in gene flow. Additionally the internal linear trunk has a large diam¬ 
eter further restricting gene flow according to the discussion of the effects 
of diameter above (compare this to the predictions of underdominance-like 
dynamics in a continuous population and the tendency of “tension zones” 
to locate in regions of restricted gene flow or population density and stop 


“pushed waves” Barton and Hewitt 

1985; [Barton and Turelli 

2011), (see 

also “invasion pinning” in 

Keitt et al. 

2001 

)• 
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3.1 Implications 

3.1.1 Natural Systems 

This work was originally motivated by biological systems. It is interesting 
to ask, based on these results, where we might expect to see underdominant 
polymorphism being maintained in wild populations. Population network 
configurations exist in which even subtle levels of underdomaince can main¬ 
tain stable geographic differences between populations with substantial rates 
of migration. Chromosomal rearrangements that lead to strong underdomi¬ 
nance can occur at relatively high rates and are rapidly established between 
closely related species (White, 1978, Jacobs, 1981). Subtle underdominant 


interactions may be more widespread than previously appreciated and may 


have played a large role in shaping gene regulatory networks (Stewart et al. 


2013). Note also that weak effects among loci can essentially self organize 


to become coupled and that this effect may extend to a broader class of 
underdominant-like effects such as the well known Dobzhansky-Muller in- 
compatabilites (Barton and De Cara, 2009; Landguth et ah, 2015). We 


have found that bifurcating tree-shaped (dendritic) networks have very high 
stability. Key natural occurrences of dendritic habitats include freshwater 
drainage systems, as well as oceanic and terrestrial ridge systems. Indeed, a 
role of underdominance in shaping patterns of population divergence across 
connected habitats has been implicated, either directly or indirectly, in fresh- 


water fish species 

Fernandes-Matioli and Almeida- Toledo[ 2001 

Alves et ah 

2003 

Nolte et ah 

2009), salamanders (Fitzpatrick et ah[[2009 

Feist, 2013), 

frogs 
rats ( 

(Bush et ah 

1977) (see also (Wilson et ah, 1974)), semiac 

uatic marsh 
Lwai‘i, which 

Nachman and Myers, 1989), and Telmatogeton flies in Hs 


rely on freshwater streams for breeding environments (Newman 


m 


the case of Dipterans we are ignoring chromosomal inversions which do not 
result in underdominance in this group Coyne et ah, 1991, 1993). In fact. 


alpine valleys around streams also follow a connected treelike branching 
pattern and there are examples of extensive underdominance in small mam¬ 
mals found in valleys in mountainous regions (Pialek et ah, 2001; Bassetj 


et ah, 2006; Britton-Davidian et ah, 2000). To the extent that persisting 


underdominant and underdominant-like fitness effects may promote specia- 


tion (rates of karyotype evolution and speciation are correlated Bush et ah 


1977) it should be noted that freshwater streams contain 40% of all fish 


species yet are only 1% of the available fish habitat and that a higher rate 


of speciation is indeed inferred for freshwater versus marine systems (Bloom 


et ah, 2013). 


However, there are also examples of the maintenance of underdominant 
polymorphisms that are not found in species associated with limnological 
structures. For example, the island of Sulawesi itself has an unusual branch¬ 
ing shape and a large number of terrestrial mammal species with a 90% or 
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greater rate of endemism excluding bats (Groves, 2001). Other factors that 


are associated with underdominant stability are the diameter of the network 
and having an even order of nodes. The Hawaiian islands essentially form 
a linear network of four major island groups (Ni‘ihau & Kaua‘i—0‘ahu— 
Maui Nui—Hawai‘i) and are known for their high species diversity and rates 


of speciation with examples in birds (Lerner et al., 2011), spiders (Gillespie 


2004), insects (Magnacca and Price 2015), and plants (Helenurm and Gan¬ 


ders, 1985 Givnish et al., 2009). For example, the Hawaiian Drosophila cras- 


sifemur complex has maintained chromosomal rearrangements between the 


islands that are predicted to result in underdominance (Yoon et al. 1975). 


Perhaps the network topology of the Hawaiian archipelago (in addition to 
the diversity of micro-climates, environments, and ongoing inter-island colo¬ 
nizations) has contributed to the high rates of diversification found on these 
islands. 

In contrast, areas where we might expect to see less maintenance of ge¬ 
netic diversity that can contribute to boundaries to gene flow are in highly 
interconnected networks with low diameters such as, perhaps, marine broad¬ 
cast spawners with long larval survival times that are associated with the 
shallow waters around islands (i.e., the network nodes). Examples of a lack 
of speciation in such groups, distributed over areas as large as half of the 


Earth’s circumference, exist (Palumbi, 1992; Lessios et ah, 2003). 


Another type of network is one that is distributed over time rather 
than space. Underdominant interactions have been inferred in the Ameri¬ 
can bellflower Campanula americana (Galloway and Etterson, 2005). This 


species is unusual in that individuals can either be annual or biennial de¬ 
pending on the time of seed germination. Given that the majority of seeds 


are expected to germinate within a single year (Galloway, 2001), even-year 


biennials may form a somewhat distinct population from odd-year biennials 
with gene flow occurring by the subset of annual plants—forming an even or¬ 
dered network. Einally, a tantalizing combination exists in the pink salmon, 
Oncorhynchus gorbuscha, of the North Pacific. This species is both biennial 
and returns to native freshwater streams to spawn. Indeed, artificial crosses 
between even and odd year individuals (of the same species) have revealed 
extensive genetic differences with hybrid disgenesis ([Gharrett and Smoker 


1991; Limborg et ah, 2014). 


3.1.2 Applications 

Various “transgene mitigation” methods have been proposed to prevent the 
transfer of genetic modifications from genetically engineered crops to tradi¬ 


et al. 


tional varieties or wild relatives (Lee and Natesan 2006; Daniell, 2002 Hills 


2007; Kwit et al., 2011) including the use of underdominant constructs 
(Reeves and Reed, 2oi4[ Soboleva et al. 2003). In a species with limited 
pollen dispersal, it may be tempting to plant a buffer crop area between 
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a GMO crop with underdominant transgene mitigation and an adjacent 
unmodified population. However, these results suggest that, over multiple 
generations, this configuration may actually destabilize the system and pro¬ 
mote the spread or loss of the genetic modification. In a simplistic scenario 
a single flanking buffer population results in an odd number of populations, 
breaking the evenness rule of stability (unless the populations are arranged 
in essentially aV = 5 star pattern, Figure [^. Depending on local condi¬ 
tions, it may be preferable to plant two distinct yet adjacent buffer crops or 
none at all. Of particular note are genetically modified and wild species that 
exist in freshwater systems, such as rice (Lu and Yang, 2009) and fish (De-| 


vlin et ah, 2006). Our predictions suggest that underdominant containment 
may in general have enhanced stability in these situations. However, this 
can work both ways. Genetic modifications may be amiable to underdomi¬ 
nant mitigation strategies to prevent establishment in the wild; yet, difficult 
to remove from freshwater systems if established. 

Using underdominance to stably yet reversibly genetically modify a wild 
population is one goal within the held of genetic pest management. The 
potential implications of these results depend on the amount of modihed in¬ 
dividuals that could be released into the wild. If the numbers are sufficiently 
high to transform an entire region then highly interconnected populations 
would be ideal to ensure full transformation. However, as is much more 
likely, if the number of individuals that can be released is much smaller 
than the total wild population, transformation might best be achieved in 
a stepwise strategy utilizing linear or treelike population conhgurations. In 
the case of Hawai‘i, limiting the effects of avian malaria by modifying non¬ 
native Gulex mosquitoes has been proposed as a method to prevent further 

2002[ ) . Linear island 


extinctions of native Hawaiian forest birds (Clarke 


archipelagos and their river valleys (Gulex are more common at lower eleva¬ 
tions, van Riper HI et ah, 1986) may be ideal cases to both transform local 
populations yet prevent genetic modifications from becoming established 
outside of the intended area. 


3.2 Future Directions 

It will come as no surprise to an evolutionary biologist that systems with 
greater genetic isolation (such as freshwater streams versus marine envi¬ 
ronments) will lead to increased genetic divergence and rates of speciation; 
however, the implication we are focusing on here is the influence of the pop¬ 
ulation network topology. We are suggesting that, for the same degree of 
migration rate isolation, alternative network topologies might be compared 
to inferred rates of speciation and/or enhanced genetic diversity that leads 
to hybrid dysgenesis. A consideration of the geological history is also ap¬ 
propriate to incorporate effects such as stream capture and the merging of 
islands on biological diversity. A proper meta-analysis or experimental evo- 
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lution of this network topology effect is beyond the scope of the current 
manuscript but would be useful projects to further explore these effects. 
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